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Graphene stacked in a Bernal configuration (60° relative rotations between sheets) differs electron- 
ically from isolated graphene due to the broken symmetry introduced by interlayer bonds forming 
between only one of the two graphene unit cell atoms. A variety of experiments have shown that 
non-Bernal rotations restore this broken symmetry; consequently, these stacking varieties have been 
the subject of intensive theoretical interest. Most theories predict substantial changes in the band 
structure ranging from the development of a Van Hove singularity and an angle dependent electron 
localization that causes the Fermi velocity to go to zero as the relative rotation angle between sheets 
goes to zero. In this work we show by direct measurement that non-Bernal rotations preserve the 
graphene symmetry with only a small perturbation due to weak effective interlayer coupling. We 
detect neither a Van Hove singularity nor any significant change in the Fermi velocity. These re- 
sults suggest significant problems in our current theoretical understanding of the origins of the band 
structure of this material. 

PACS numbers: 73.22.Pr, 61.48.Gfi, 79.60.-i 
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I. INTRODUCTION 

Multilayer epitaxial graphene (MEG) grown on the 
SiC(OOOi) (C-face) is now known to have a highly-ordered 
non-Bernal (non-AB) stacking, where adjacent graphene 
planes have commensurate relative rotations that are not 
60°iii^ These commensurate rotations lead to large su- 
percells [as shown in Fig. [Ija)] that are seen in STM 
images of MEG films [see Fig. ^h)]. These STM viewed 
supercells are often referred to as "Moire" patterns. 

Since the discovery of C-face MEG films, the effect 
of this new graphene stacking on the band structure of 
these films has been the subject of active experimental 
and theoretical study. The earliest ab-initio calculations 



for a large-angle commensurate graphene bilayer rotation 
(32.20°) predicted essentially no effect on the graphene 
band structure near the _fC-points of either graphene 
sheetjii^ In other words, the band structure consists of 
two independent but rotated Brillouin Zones (BZ) as 
shown in Fig. [IJc) where nearby Dirac cones from the 
two sheets at K and Kg do not interact as shown in 
Fig. I^a). This prediction has been borne out in a num- 
ber of experiments on MEG films. Electron transport}^ 
infrared adsorption spectroscopy,^ angle-resolved pho- 
toemmision (ARPES))^'^ and scanning tunneling spec- 
troscopy (STS)'* all show that the graphene sheets in 
these films behave nearly identically to electronically- 
isolated graphene sheets. 
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FIG. 1. (a) Schematic of the commensurate rotation of two graphene sheets. The commensurate angle 6 leads to a large 
supercell structure, (b) a 200 x 200A STM image of a {m,n) = (4,5) supercell for C-face graphene sheets with a relative 
rotation of 6 = 7.34° ^ (c) The hexagonal Brillouin Zone of two graphene sheets with a relative rotation 6. The smaller Brillouin 
Zone of the commensurate superlattice (defined by reciprocal lattice vectors Gi and G2) is also shown. 




FIG. 2. (a) Schematic Dirac cones of two commensurately 
rotated graphene sheets, The cones are shifted by a vector 
Kg. (b) Velocity renormalization causes vp to be smaller than 
an isolated sheet making the slopes of the cones smaller in an 
energy window AEj^ near the Dirac point, (c) Predicted Van 
Hove singularity formed between the K- and /i'e-points of the 
two rotated Dirac cones. 



Numerous theoretical treatments for rotated graphene 
bilayers, however, have predicted more dramatic effects 
in the band structure of these films. Continuum ap- 
proximations to a tight-binding model, as well as full 
tight-binding (TB) and ab inito calculations on small an- 
gle commensurate graphene bilayers, predict a substan- 
tial decrease (renormalization) in the Fermi velocity, vp, 
near the ii'-point (Dirac point). In these calculations the 
interlayer interaction causes the wave functions in the 
graphene sheet to become highly localized for small rela- 
tive rotations of the bilayeri^i These localized states have 
a reduced vp that is exceptionally small when the rela- 
tive rotation is less than 3° - less than 25% of the value 
for an isolated graphene sheet. This velocity reduction 
occurs in an energy window, ^Efr, defined as the energy 
difference between the Dirac point energy (i.e. at the K- 
point) and the energy where the two rotated Dirac cones 
cross (the Af '-point) [see Figs.IlJc) and HJb)] J^^i^ 



In addition to the velocity renormalization. Lopes dos 
Santos, et alJ^ have also predicted that the interlayer 
interaction leads to an angle dependent Van Hove singu- 
larity in rotated bilayers. The singularity appears where 
the Dirac cones from the two rotated layers cross at the 
Af '-point [see Figs.IlJc) andHJc)]. 

MEG films offer a perfect platform for definitive tests 
of these predictions. While MEG films have average rel- 
ative rotations of ~30°, much smaller relative rotations 
occur with significant frequency.— i^i^ Furthermore, these 
films are flat^ extremely well-ordered and can be 
grown as thick multilayers or as thin as single sheets. 

In this work, we use high resolution ARPES to di- 
rectly measure the band structure of commensurate ro- 
tated graphene sheets. These experiments conclusively 
show, despite the theoretical predictions for large reduc- 
tions vnvp and the development of a Van Hove singular- 
ity, that commensurately-rotated graphene sheets show 
no significant deviations from the linear band structure 
of graphene. Measurements on many samples and many 
relative rotations, from both thick and thin films, show 
that the band structure of MEG films remains nearly 
identical to that of an isolated graphene sheet. 



II. EXPERIMENTAL 

The substrates used in these studies were n-doped 
n = 2 X lO^^cm"^ 6H-SiC. Samples were grown in a 
closed RF induction furnace using the Controlled Silicon 
Sublimation (CSS) method.^"* The samples were trans- 
ported in air before introduction into the UHV analysis 
chamber. Prior to ARPES measurements the graphene 
films were thermally annealed at 800 °C in UHV. Sample 
thickness was measured by ellipsometryJ^ ARPES mea- 
surements were made at the Cassiopee beamline at the 
SOLEIL synchrotron in Gif/Yvette. The high resolution 
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Cassiopee beamline is equipped with a modified Peterson 
PGM monochromator with a resolution Ej t^E ~ 70000 
at 100 eV and 25000 for lower energies. The detector is a 
±15° acceptance Scienta R4000 detector with resolution 

A£;< ImeV and Afc~O.OlA"^ at ?iw = 36 eV. AU mea- 
surements were carried out at 4K. The total measured 
instrument resolution is (A_E<12meV). 



III. RESULTS 

To understand subsequent ARPES data, we first re- 
view how commensurate rotations determine the band 
structure of rotated graphene. Using the notation of 
Melefii a commensurate rotation is determined by the 
supercell vector a = mt\ -\- nt2, where m and n are inte- 
gers and ti and are the unit cell vectors of graphene. 
The commensurate relative rotation of the two sheets is 
determined by integers ra and n; cos Q — [Aran -\- t? -V 
m^) / 2{m'^ +mn) . In reciprocal space the commensu- 
rate supercell forms a small Brillouin Zone (BZ) defined 
by the two reciprocal lattice vectors Gi and G2 defined 
by; 



^Jl — ^gl 



-Q 



92 



Tl TfX 
G2 = Qgl H Qg2, 

C C 



(la) 



(lb) 



where c = rn^-f n^-t-ron and Qgi and Qg2 are the reciprocal 
lattice vectors of the unrotated graphene. This leads to 
a supercell rotated by an angle cf) relative to Qgi; cos (f) — 
(m+nj'i) I \/rrP- -I- + ran. The reciprocal space picture 
then consists of two graphene BZs with a relative rotation 
d as shown in Fig. [TJc). 

Because the two sheets are commensurate, the K- 
points of the two BZs can be connected by linear combi- 
nations of commensurate supercell reciprocal lattice vec- 
tors, i.e. Q = pGi+qG2, where p and q are integers/^ De- 
pending on the symmetry of the supercell either K— Kg = 
AKg = g (SE-even symmetry!^) or K'-Ke = AK^, = g 
(SE-odd symmetr jiii) , but not both simultaneously [see 
Fig. [lJc)]J^ The g vectors are related to the rotation 
angle by 



\g\^\AKe\^2KrKsm9/2, (2a) 



\g\^\AK'g\^2KrKsm{6O-0)/2. (2b) 

Where Kyk = \Qg\/Vi- Without considering interac- 
tions between sheets in a bilayer, the adjacent Dirac cones 
at K and Kg are simply interleaved cones as shown in 
Fig. Ha). 

Going beyond independent graphene sheets is the sub- 
ject of many theoretical works. For small rotation angles 



Lopes dos Santos, et aliiS have made an analytic solu- 
tion to the tight-binding model for small relative rota- 
tions. They use a continuum approximation to describe 
the large range of interlayer bonding geometries in the 
large supercells that result from a relative rotation. By 
assuming that the interlayer hopping has a long wave- 
length nature so that the coupling between different val- 
leys {K and K') can be ignored, they treat the interlayer 
coupling as being uniform. This calculation predicts two 
observables. First, the two nearby cones separated by g 
re-hybridize to form a distorted single cone with a Van 
Hove singularity (saddle point) halfway between the two 
K and i^e-points as shown in Fig. [2Kc). The second pre- 
diction is that the Fermi velocity slows to a renormalized 
value, vp^ that scales with rotation angle approximately 



as 



10 



(3) 



where t±^ is the interlayer hopping parameter for the ro- 
tated bilayer. For small rotations angles is assumed to 
be independent of angle (<_L~0.4i_L, 0.3eV),i^ While 
Eq. [3] overestimates the reduction mvp below ~2°, both 
full TB and ab initio calculations confirm the continuum 
approximation prediction that the renormalized velocity 
goes to zero very fast at rotation angles less than 4° 

Both the existence of a VHS and a renormalized vp 
can be tested by directly measuring the band structure of 
MEG films. To make such comparisons, the band struc- 
ture from rotated commensurate graphene sheets with 
known relative rotations must first be identified. We now 
describe exactly how commensurate rotation angles are 
measured in ARPES. 

Within the detector area of an ARPES measurement 
at huj = 36eV, Dirac cones from graphene sheets with rel- 
ative rotations angles of ±15° are visible. Some of the 
visible cones come from different uncorrelated rotated do- 
mains, i.e from a continuous graphene film that has an 
induced rotational boundary, A<j), caused by pleats in the 
graphene or by substrate step edges, etc. that are within 
the 30/im beam diameter. To distinguish the commen- 
surate pairs formed by stacked rotated graphene sheets 
from incoherent pairs formed from pleats and steps, etc., 
we take advantage of the fact that in ARPES the photo- 
emitted electron can diffract from the local surface struc- 
ture. If there is no superlattice, as in the case of cones 
from two incoherent rotational domains, there will be 
no diffraction. On the other hand if the cones are from 
a coherent bilayer pair that form a superlattice, addi- 
tional replica cones caused by diffraction will be visi- 
ble in ARPES. This diffraction effect has been nicely 
demonstrated for graphene grown on the Si-face of SiC 
where a (6-\/3 x 6-\/3)R30° superstructure causes replica 
Dirac cones at ^-vectors relative to the main graphene 
Dirac cones [see Refs. [l^ and[2d|. In the case of C-face 
graphene, the APRES signature of two stacked commen- 
surately rotated graphene sheets will not only be the two 
primary Dirac cones but will also include replica (diffrac- 
tion) cones positioned at ±g relative to either of the pri- 
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FIG. 3. (a) The band structure near the K- and _K'e-points {kx — fcrx) of two closely rotated planes (marked by arrows), ky is 
perpendicular to the F — K direction. The relative rotation angle between the two cones is 4.2°. (b) a schematic of the BZ of 
two rotated graphene sheets with a blowup near the K and Kg points. Two different detector planes (dashed arcs) are shown 
cutting through the Dirac cones and their replicas at itnQ. Note that the constant energy cuts through the cones are not full 
circles as described in text, (c) a similar ARPES scan of (a) but taken at a detector cut slightly closer to the F point. Data is 
plotted on a log scale to show the supercell C/-replicas band (R) (marked by a third arrow). 



mary cones [ see Fig.[3][b)]. 

An example of how a commensurate pair is identified is 
shown in Fig. [3] FigurelSja) shows a MEG band structure 
measured by ARPES near the graphene K- and Kg-point 
of two closely spaced Dirac cones (A/cj^ = 0.124 ± .OOSA, 
(kx — krK ) ■ The two cones are identified as being part 
of a commensurate rotated pair as demonstrated in the 
log intensity plot of Fig. ^c). The data is plotted in 
a log scale because the diffraction cross-section of the 
photo-electron is weak. This makes the replica intensity 
significantly smaller compared to the primary cones 
The kx distance to the F point in Fig. [2^c) is slightly 
reduced relative to the image in Fig. ^a) to enhance the 
replica's intensity [see Fig. ^h) ]. The replica cone R 
(marked by a third arrow in Fig. [3fc)) has the same sep- 
aration from the K point as the K point is from the Kq 
point. Note that because the detector plane is not exactly 
perpendicular to the TK direction and because constant 
energy cuts through the cones are only partial circleSj^ 
replicas appear more intense in the +ky direction. Using 
Eq. [5J the C/-vector connecting the pair of Dirac cones in 
Fig. [3ja) allows us to find the relative rotation of the two 
sheets to be 4.15 ± 0.07°. Note that the fc^-width of the 

bands is resolution limited (~ O.OlA ^), indicating that 
the supercell domain size is well ordered over a length 
scale of 800A. 

We can be more restrictive in how we identify com- 
mensurate rotations. In a commensurate pair, one of the 
graphene pairs must lie below the other. At the photon 
energy used in these experiments, the electron mean free 
path. A, is only 3.4A (approximately the graphene in- 
terlayer thickness). This means that one cone in a pair 
should have relative intensity that is ~ 40% of the most 



intense cone in the pair. For the analysis used in the 
work, we analyze only commensurate rotated pairs iden- 
tified by both the existence of a related replica cone and 
that the two cones have the proper relative intensity for 
bilayers. 

We have used ARPES spectra from similar cones pairs 
as that shown in Fig. [3] to investigate predictions for ve- 
locity renormalization vp- The Fermi velocity is derived 
from the slope of Dirac cones for relative rotations as 
low as 1.1°, for multiple samples, and for film thicknesses 
ranging from 3-10 layers. All velocities were measured in 
an energy window between E]j and the energy where the 
pair of cones cross (the M'-point in Fig. [T] . The results 
are plotted in Fig. [H 

It is clear from Fig. H] that no significant velocity 
changes are observed even for the smallest measured ro- 
tation angles, where the velocity has been predicted to be 
50 times smaller than the velocity for isolated graphene. 
It must be pointed out that no cones, regardless of wether 
or not replica cones were visible, had any detectable ve- 
locity renormalization. The average (pp) for all mea- 
sured cones from several samples and graphene thick- 
nesses is 0.99 ±0.05 X lO^m/sec. Within error bars, the 
measured (wp) is consistent with values obtained for exfo- 
liated graphene^^ and from both IR measurementSi^ and 
STS^ on MEG films. 

While the predicted velocity renormalization is a prob- 
lem by itself, we are also able to show that no Van Hove 
singularity forms at the M'-point between the two cones 
[shown schematically in Fig. [5Jc)]. Figure [5] shows two 
detailed view near the crossing of two commensurately- 
rotated cones. The cones cross as straight lines with 
no singularity at the crossing (i.e, the Af '-point). The 
expected TB dispersion with VHS is drawn in Fig. [5] 
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Relative Rotation 8 

FIG. 4. Comparison of the measured vp (o) versus the 
predicted renormalized velocity vf as a function of rela- 
tive rotation angle. Solid line is the approximate vf from 
Ref. [13]. A's are the low angle corrected vf from a full TB 
calculation.— Dashed line is the average velocity (wf) for all 
measured Dirac cones. 

for clarity. To within the experimental energy and mo- 
mentum resolution, we can say that the crossing is an 
undistorted intersection of two linear bands. Numerous 
samples and cones have been measured and no VHS are 
observed even for small angles like those in Fig. EJb); 
precisely where the continuum model should be most ap- 
propriate. Furthermore, no increase is observed in the 
integrated spectral density (proportional to the density 
of states) at the crossing that would indicate a change 
in \/i^{E) [this is discussed in more detail in Sec. llVj . 
The experimental uncertainty places an upper limit on 
the size of a possible gap at the crossing to be less than 
60meV. While this is well below predicted values for a 
the Van Hove singularity of Lopes dos Santos, et aL^^ 
the uncertainty does not rule out the small gaps pre- 
dicted to occur for a class of rotations with weak inter- 
layer couplingii^ 



IV. DISCUSION 

The ARPES measurements presented in Sec. lIIII clearlv 
demonstrate serious inconsistencies between the calcu- 
lated and experimental band structure from commensu- 
rately rotated graphene sheets. Both the predicted Van 
Hove singularity and renormalized Fermi velocity are not 
observed, even for small relative rotations. While we 
cannot offer definitive reasons for the differences between 
current theory and these experiments, we can provide ad- 
ditional data that places limits on any model that may 
be used used to reconcile these differences. 

First, it must be emphasized that surface effects can- 
not be the source of the inconsistency between theoreti- 
cal predictions for vp and the values measured in these 
experiments. This is because, even in thick MEG films 
(fO— 40 layers), infrared absorption experiments measure 
the same velocity;^ 

A more important consideration is the valid energy 
range where renormalization is significant. Models dif- 
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FIG. 5. A close up view of the band structure from (a) 4.2° 
and (b) 1.7° commensurate graphene pairs. The dashed line 
is the predicted band structure with a Van Hove singularity 
from a continuum TB model.— Note that the cones are hole 
doped by ~ 15meV. 

fer in the energy range where the renormalized velocity 
prediction is expected to be valid. More rigorous cal- 
culations, beyond first order expansions to a TB model, 
predict that at A:-values larger than the first supercell BZ, 
the velocity should return to the isolated graphene value 
- assuming no Van Hove singularity forms at the M'- 
point.i^ To ensure that the experimentally measured ve- 
locities are analyzed in an energy region valid for all mod- 
els, we need to choose a minimum energy below the Dirac 
point where measured velocities are insured to be valid. 
The predicted renormalized vp should be exact along the 
TK direction within the first BZ of the commensurate 
supercelliiiii^ The energy window below the Dirac point 
is then set by AKe and vp; AE^k = ^vp\A'Kg\. A more 
restrictive window would be to use the energy of the band 
at the cone crossing M'-point [see Fig. [Tfc) ], AEkm- 
Using Eq. [31 the energy where the renormalized band 
intersects the M'-point would be: 

AEkm « " 9{i±/nvpo\AKg\n (4) 

We plot the normalized slope of the Dirac cones, 
{l/h)[dE/dk){\/vpo) = Vp/vpo, in Fig. [Hi as a function 
of energy for rotated graphene pairs for three different 
rotation angles. The energy AEkm and the theoretical 
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E-E^ (eV) 

FIG. 6. Measured velocity, derived from the slope of the Dirac 
cones, as a function of energy. The relative rotation angles 
are (a) ^ = 4.1° and (b) ^ = 2.4°. The velocity is normalized 
by vfo ~ lO^m/sec. The horizontal dashed line marks the 
TB prediction from Eq. [31 Vertical line marks the predicted 
renormalized band energy AEkm at the supercell M-point. 
The shaded region is the energy regime where the two cones 
cross and where E{k) is subsequently difficult to accurately 
determine. 



velocity reduction are marked for each rotation angle in 
Fig. [5] For the three angles shown, the measured veloc- 
ity between AEkm and Ep is always higher than the 
predicted TB value and nearly equal to lO^m/sec. We 
note that the change in slope close to the Dirac point 
('^ lOmeV) should not be over interpreted. The energy 
resolution coupled with the fc-resolution (0.01 A ) and 
the distortion due to the Fermi-Dirac function near Ep 
mean that it is difficult to accurately fit the position of 
the two converging tt bands. A more important observa- 
tion is that the measured velocity does not go to zero at 
AEkm as it should if there was a Van Hove singularity 
at the M'-point. These results show that if any signifi- 
cant velocity renormalization occurs, it is on an energy 
scale much smaller than current predictions. 

In contrast to this work, a singularity in the dl/dV 



curves of STS experiments from CVD grown graphene 
deposited on graphite has been interpreted as evidence 
of a Van Hove singularity caused by rotated graphene 
sheets, although no corresponding velocity change has 
been reported. ^'^ It is worth noting that there is a signif- 
icant difference between MEG and graphene on graphite 
that makes the applicability of all the theoretical models 
discussed above tenuous for the graphene on graphite sys- 
tem. First, the Moire patterns of graphene on graphite 
are not necessarily due to simply sheet rotations. For ex- 
ample, surface graphite Moire patterns are the result of 
a mechanically distorted top layer formed during cleav- 
ing of graphite. Similar distortions and buckling are 
seen in depositing exfoliated graphene onto a support 
substrate. The consensus understanding is that these 
Moire patterns are not caused by the relative rotation of a 
graphene sheets>24 In other words, it is possible to have a 
Moire pattern without rotations. Furthermore, the Moire 
corrugation amplitude in graphene on graphite films is 
typically 1 — 2 A [see Pong and Durkan^i and references 
therein]. This should be compared to MEG films where 
the corrugation is less than 0.2A,^-^'^ Height modulations 
that are more than a third of the graphene interlayer 
are certainly not part of the any of the models discussed 
above and therefore make a comparison between theories 
based on flat rotated sheets questionable. 

It has also been argued that the lack of a Van Hove 
singularity in MEG films implies that there is no inter- 
layer coupling {t± ^0).^4 This conjecture oversimplifies 
the problem for a number of reasons. First, the cou- 
pling must be non-zero otherwise the films would simply 
delaminate in solution, which they do not. Second, it 
is theoretically possible in TB to have a coupling and 
no Van Hove singularity. Mele has recently shown that 
a certain class of rotation angles (those with SE-even 
rotationaii ) , while having a non-zero interaction, do not 
have a singularity at the supercell M'-point.— Third, re- 
cent STS experiments on the fine structure splitting of 
the zero-Landau level in MEG films indicate that third 
layer interaction are consistent with TB interlayer cou- 
pling parametersi^ 



V. CONCLUSION 

We have shown by directly measuring the band struc- 
ture of multilayer epitaxial graphene that serious incon- 
sistencies exist between the calculated and experimental 
band structure from commensurately rotated graphene 
sheets. Both the predicted Van Hove singularity and the 
expected large reduction in the Fermi velocity, as the 
commensurate rotation angle between graphene sheets 
goes to zero, are not observed, at least not at the dra- 
matic level predicted or in the energy window specified 
by current theories. Given that these theoretical predic- 
tions are broad-based, encompassing both tight-binding 
and ab intio calculations, the inconsistency between the- 
ory and experiment suggests a fundamental process that 
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is not included in current theoretical treatments. It is 
possible that local strain fields caused by nearly overlap- 
ping TT-bonds in the commensurate supercell lead to small 
relaxations that break some of the symmetry of the ideal 
rotated graphene pair. While such speculations are inter- 
esting, they lie outside the scope of this work. Nonethe- 
less, the discrepancies pointed out this work suggest that 
a detailed look at these problems will be an important 
theoretical research avenue. 
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